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We present detailed analysis of the two gamma-ray sources, IFGL J1801.3-2322c 
and IFGL J1800.5-2359c, that have been found toward the supernova remnant (SNR) 
W28 with the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space 
Telescope. IFGL J1801.3-2322c is found to be an extended source within the 
boundary of SNR W28, and to extensively overlap with the TeV gamma-ray source 
HESS J1801-233, which is associated with a dense molecular cloud interacting 
with the supernova remnant. The gamma-ray spectrum measured with LAT from 
0.2-100 GeV can be described by a broken power-law function with a break of 
~ 1 GeV, and photon indices of 2.09 ± 0.08 (stat) ± 0.28 (sys) below the break and 
2.74 ± 0.06 (stat) ± 0.09 (sys) above the break. Given the clear association between 
HESS Jl 801-233 and the shocked molecular cloud and a smoothly connected spec- 
trum in the GeV-TeV band, we consider the origin of the gamma-ray emission in 
both GeV and TeV ranges to be the interaction between particles accelerated in the 
SNR and the molecular cloud. The decay of neutral pions produced in interactions 
between accelerated hadrons and dense molecular gas provide a reasonable explana- 
tion for the broadband gamma-ray spectrum. IFGL J1800.5-2359c located outside 
the southern boundary of SNR W28, cannot be resolved. An upper limit on the size 
of the gamma-ray emission was estimated to be ~ 16' using events above ~ 2 GeV 
under the assumption of a circular shape with uniform surface brightness. It appears 
to coincide with the TeV source HESS J1800-240B, which is considered to be associ- 
ated with a dense molecular cloud that contains the ultra compact H II region W28A2 
(G5. 89-0. 39). We found no significant gamma-ray emission in the LAT energy band 
at the positions of TeV sources HESS J1800-230A and HESS J1800-230C. The LAT 
data for HESS J1800-230A combined with the TeV data points indicate a spectral 
break between 10 GeV and 100 GeV. 

Subject headings: cosmic rays — acceleration of particles — ISM: individual objects 
(W28, G6. 4-0.1) — ISM: supernova remnants — gamma rays: ISM 



1. Introduction 

Diffusive shock acceleration operating at supernova shock waves can distribute particles to 



very h igh energies with a power-law form having number index about 2 (e.g. jBlandford & Eichlei 



19871). It is generally expected that if a dense molecular cloud is overtaken by a supernova blast 



wave, the sho cked molecular cloud c an be illuminated by relativistic particles accelerated at super- 



nova shocks (|Aharonian et al.lll994|) . If the accelerated particles are comprised mostly of protons. 



say > 100 times more abundant than electrons like the observed Galactic cosmic rays, decays of 



neutral pions produced in inelastic collisions of the accelerated protons with dense gas are expected 
to be a dominant radiation component in the gamma-ray spectrum of the cosmic-ray-illuminated 
molecular cloud. Although an earlier attempt to detect TeV gamma-ray emission from supernova 
remna nts (SNRs) that have evidence for molecular cloud interactions with the Whipple telescope 
failed (jBu cklev et al.lll998b. two arch etvpical S NRs interacting with m olecular clouds, IC 443 
jAlbert et al...2007l : lAcciari et alJboog I) and W28 jAharonian et allboosl) . have been detected with 
the current generation of imaging atmospheric Cherenkov telescopes. However, the identities of 
the particles responsible for the TeV sources remains elusive. 

The advent of the Large Area Telescope (LAT) onboard the Fermi Gamma-ray Space Tele- 
scope has brought a new opportunity to study the gamma-ray emission from SNRs at G eV energies. 



LAT observa tions of middle-ag ed SNRs inte racting with molec ular clouds, W51C (|Abdo et al 



20091), W44 (lAbdo et al.ll2010ah . and IC 443 (lAbdoet all 1201 Obh . have revealed bright extended 
gamma-ray sources coincident with the SNRs. The gamma-ray luminosity reaches ~ 10^^ erg s~\ 
which immediately rules out an inverse-Compton origin of the GeV gamma rays since it requires 
total electron energy comparable to or larger than the typical kinetic energy released by a supernova 
explosion, ~ 10^^ erg. The gamma-ray spectra of the three remnants exhibit remarkable spectral 
breaks at an energy of several GeV, making these SNRs much less luminous at TeV energies. This 
characteristic demonstrates the importance of observations at GeV energies. 

W28 is a mixed-morphology SNR, characterized by center-filled thermal X-ray emission and 
shell-like radio morpholo gy. In addition X-ray ob servations show limb-brightened shells in the 
northeast and southwest (IRho & Borkowskil |2002|) . The shell-like radio emission is prominent 



in the northeastern region with slightly fainter emission a t the northern b oundaries (|Dubner et al. 



2000|) . Interactions of the SNR with molecular clouds dWoottenlll98lh along its no rthern anc 
northeastern boundaries are traced by the high concentration of 1720 MHz OH mas ers (IFrail et al. 



1994: Claussen et al. 



1997l . ll999h . and high density (hh > 10^ cm ^) shocked gas (lArikawa et al. 



19991 : lReachetal.ll2005b . The overall shape of W28 is elliptical with a size of 50' x 45'. W28 
is located within a complex, star-forming region along the Gal actic plane to ward the large H II 
regions (M8 and M20) and young clusters (e.g., NGC 6530) (iGoudisI 19761). The o bservations 



of molecular lines place SNR W28 at a distanc e of ~ 2 kpc (I Velazquez et al. Il2002() . Estimates 
for its age vary between 3.5 and 15 x 10"* yrs (IKaspi et al.lll993h. W28 is considered to be an 



evolved remnant in th e radiative stage o f SNR evolution (|Lozinskay all 1 992() . which is supported by 
optical observations dLozinskay al 1 1 974|) . Measurements with Energetic Gamma-Ray Experiment 
Telescope (EGRE T) onboard the Comp ton Gamma-ray Observatory found a gamma-ray source, 
3EG J 1800-233 8 faartman et al.lll999l) in the W28 field. However, its association to SNR W28 
remained unclear mainly due to large source location uncertainties from EGRET. A gamma-ray 
source is listed i n the W28 field in the AGILE (Astro-rivelatore Gamma a Immagini LEggero) 
one year catalog (|Pittori et al.ll2009|) . However, detailed analysis of this field is not published by 
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AGILE yet. 

H.E.S.S. observations of the W2 8 field have revealed fo ur TeV gamma-ray sources position- 
ally coincident with molecular clouds (lAharonian et alJlOOSh : HESS Jl 801-233, located along the 



northeastern boundary of W28, and a complex of sources, HESS J1800-240A, B and C, located 
~ 30' south of SNR W28. HESS J1801-233 coincides with a molecular cloud interacting with 
SNR W28, providing one of the best examples of a cosmic -ray-illuminated cloud. Understanding 
the origins of TeV emission in HESS J 1800-240 ABC is of particular interest; they may be due to 
cosmic rays that have diffused from W28. 

In this paper, we report Fermi LAT observations of gamma-ray sources in the W28 field in 
the GeV domain. First, we give a brief description of the observation and gamma-ray selection 
in Section [2l The analysis procedures and results are explained in Section [3l where the spatial 
extension and spectra of the LAT sources in the W28 field are described. Discussion is given in 
Section in followed by conclusions in Section [51 



2. OBSERVATION AND DATA REDUCTION 

The LAT is the main instrument on Fermi sensitive to gamma rays. The energy range spans 
from ~ 20 MeV to > 300 GeV, although as noted below in the present analysis we use only events 
with energies > 200 MeV. It is an electron-positron pair production telescope, built with tungsten 
foils and silicon microstrip detectors to measure the arrival directions of incoming gamma rays, and 
a hodoscopic cesium iodide calorimeter to determine the gamma-ray energies. They are surrounded 
by 89 segmented plastic scintillators that serve as an anticoincidence detector to reject charged 
particle eve nts. Details of the LAT instrument and pre-launch expectations of the performance can 



be found in lAtwood et al.l (120091) . Relative to earlier gamma-ray missions, the LAT has a large 
~ 2.4 sr field of view, a large effective area (~ 8000 cm^ for >1 GeV if on-axis) and improved 
angular resolution or point-spread function (PSF, better than 1° for 68% containment at 1 GeV). 

Routine science operations of the LAT began on August 4, 2008, after the conclusion of a 
commissioning period. We have analyzed events in the W28 field, collected from August 4, 2008, 
to July 4, 2009, with a total exposure of ~ 2.8 x 10^° cm^ s (at 1 GeV). During this time interval, 
the LAT was operated in sky survey mode nearly all of the time. In this observing mode the LAT 
scans the sky, obtaining complete sky coverage every 2 orbits (~ 3 hr) and relatively uniform 
exposures over time. 

We used the standard LAT analysis software, ScienceTools version v9rl5, which is available 
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from the Fermi Science Support Center (FSSCU and applied the following event selection criteria: 
a) events have the highest probability of being gamma rays, i.e., they should be classified as so- 



called Pass 6 diffuse class (JAtwood et al.ll2009|) . b) the reconstructed zenith angles of the arrival 



direction of gamma rays should be less than 105°, to minimize contamination from Earth albedo 
gamma rays, c) the center of the LAT field of view should be within 39° from the zenith in order 
to exclude data from the short time intervals when the field of view can be partly occulted by 
the earth. There are no gamma-ray bursts detected by the LAT within 15° of the W28 field, thus 
we did not need to apply any additional time cut. The energy range analyzed here is restricted 
to > 200 MeV to avoid possible large systematic uncertainties due to the strong Galactic diffuse 
emission near the Galactic center, smaller effective area, and much broader PSF at lower energies. 



3. ANALYSIS AND RESULTS 

3.1. Source position and Spatial Extension 

Figure [H shows a smoothed count map in the 2-10 GeV energy band in a 10° x 10° region 
around W28. Figure [2]gives comparisons with images of other wavebands in close-up view. Black 
contours indicate the H.E.S.S. significance map in (a), CO (J=l-0) line intensity taken by NAN- 
TEN for the velocity range from km s"^ to 20 km s "' (corresponding to kinematic distances of 
approximately to 4 kpc) in (b), and a VLA image in (c). Correlations between GeV gamma rays 
observed by LAT and some of the H.E.S.S. sources are evident. Gamma rays are also bright in the 
brightest spots in the CO contours, and bright regions in the gamma-ray image extensively overlap 
with bright reg ions in the VLA co ntours. There are two LAT sources in the vicinity of W28 in the 



IFGL catalog (lAbdo et al.ll2010dl) : IFGL J1801.3-2322c, and IFGL J1800.5-2359c. No obvious 
pulsations of gamma rays are found in these sources. Hereafter we refer to IFGL J1801.3-2322c 
as Source N, and to IFGL J1800.5-2359c as Source S. 

In order to quantitatively evaluate the extension and location of these two sources, we apply 
the maximum likelihood tool, gt like, which is publicly available as part of Fermi ScienceTools. 
The likelihood is the product of the probability of observing the gamma-ray counts of each spatial 
and energy bin given the emission model, and the best parameter values are estimated by maximiz- 
ing the likelihood of the data given the model (IMattox et al.ll996b . The probability density function 



for the likelihood analysis included a) individual sources detected in the preliminary LAT 1-year 
catalog, b) the Galactic diffuse emission resulting from cosmic-ray interactions with the interstellar 



'Software and documentation of the Fermi ScienceTools are distributed by Fermi Science Support Center at 
[http://fermi.gsfc.nasa.gov/ssc 
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medium and radiation based on the LAT standard diffuse background model gll_iem_v02 available 
from FSSCq and c) the isotropic component to represent extragalactic and residual cosmic -ray 
backgrounds using the isotropic spectral template isotropic_iem_v02 from FSSC. Note that we 
make energy dependent corrections of the the Galactic diffuse model by multiplying a power law 
function with the spectral index free to vary in the fit. This correction gives better spectral fits by 
taking account of local systematic discrepancies between the data and the Galactic diffuse model. 
The region of interest for the binned maximum likelihood analysis was a square region of 20° x 20° 
centered on W28 with a pixel size of 0.°1. The instrument response functions (IRFs) used in 
our work were the "Pass 6 v3" (P6_V3) IRFs, which were developed following launch to address 
gamma-ray detection inefficiencies that are correlated with background rates. 

Using the tools described above, we investigated spatial extensions of the two LAT W28 
sources. Here we used only events above 2.15 GeV to take advantage of narrower PSF in the higher 
energy band. Before investigating the detailed extension, we first determined the flux and spectral 
index for all components except for Source N and Source S. In this process, the normalization of 
the Galactic diffuse emission and the flux and spectral index of power-law model for the sources 
within 5 degrees of Source N were set free to account for the effects of sources around W28 
on the background flux in the fit. The flux and spectral index except for those of Source N and 
Source S are fixed hereafter. We used a radially symmetric uniform disk to evaluate the extension. 
We varied the radius and location of the disk while holding the position of Source S fixed at the 
catalog position, and evaluated the resulting maximum likelihood value (Lex) with respect to the 
maximum likelihood for no source hypothesis (Lq) and the point source hypothesis (Lps). A two- 
source hypothesis (L2s) is also evaluated by scanning positions of two point sources. 

The best likelihood ratio -21n(Lps/Lex) ~ 76 is obtained for Source N with a disk radius of 
a = 0.°39, which rejects a point source hypothesis at more than 8 a. Note that we v alidated this 



proced ure by applying this method on a nearby gamma-ray pulsar, PSR J 1809-2332 (lAbdo et al. 



2010cl) (~ 2° away from W28), where we find the extension to be consistent with a point source. 
In addition, the best likelihood ratio -21n(Lo/Lex) ~ 336 for the disk shape is substantially better 
than that for two point sources, -21n(Lo/L2s) ~ 276, where the positions of the two sources were 
free in the optimization. Therefore we use the best-fit disk model for further analysis of Source N. 
The best location of the disk model is found at (R.A., Dec.) = (18''0r"2P, -23°26'26") with an 
error radius of 0.°03 at 68 % confidence level. 

The extension of Source S is also investigated using the same procedure as above. We did not 
find significant extension. An upper limit on the size of the gamma-ray emission was obtained by 
investigating the decrease of the likelihood with increasing source size. Under the assumption of a 



^The model can be downloaded from|http://fermi.gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.litniil 



disk shape, the upper limit amounts to 16' at the 68% confidence level. We treat Source S as a point 
source hereafter. The best-fit location of Source S was estimated to be (R.A., Dec.) = (18''00'"59\ 
-24°11'31") using gamma rays above 2.15 GeV. This location is 13' away from the location of 
IFGL J1800.5-2359c using gamma rays above 0.1 GeV and assuming IFGL J1801.3-2322c is 
a point source. In order to investigate the origin of the position difference, we performed local- 
ization with the same energy cut as the catalog analysis assuming a point source hypothesis and 
the extended source hypothesis. We found that the localization with a point source hypothesis is 
consistent with the catalog position while the extended source hypothesis yields the location con- 
sistent with our localization with £ > 2.15 GeV cut. From this comparison, we conclude that the 
localization difference is mostly due to the extension of Source N. We use the average of two lo- 
calizations, (R.A., Dec.) = (18''00'"47'*, -24°05'45"), for the following analyses, with a systematic 
uncertainty of half of the positional difference, ~ 6'. 

We find no strong gamma-ray emission in the LAT energy band at the locations of HESS J 1 800-240 A 
and C. We evaluated the likelihood of point sources at these locations and found likelihood ratios 
-21n(Lo/Lps) ^ 10 (corresponding to a ~ 3 a) for HESS J1800-240A and -21n(Lo/Lps) ^ 9 (cor- 
responding to a ~ 3 cr) for HESS J1800-240C, which confirms no significant emission at these 
locations. 



3.2. Energy Spectrum 

We used the maximum likelihood fit tool, gtlike, for the spectral analysis of LAT sources 
in this region. Flux points were obtained by performing the maximum likelihood analysis in each 
energy bin. We used eight logarithmically spaced energy bins from 215 MeV to 100 GeV for 
Source N. We did not divide the energy range from 0.2 GeV-100 GeV into bins for Source S due to 
poor statistics. Fig.[3]shows resulting SEDs for (a) Source N and (b) Source S. The 68% confidence 
region is illustrated for Source S up to ~ 20 GeV. We do not have sufficient number of events above 
~ 20 GeV to constrain the spectrum. Upper limits at 90% confidence level are calculated assuming 
a photon index of 2 if the detection is not significant in an energy bin, i.e., the likelihood ratio with 
respect to no source is less than 9. Note that the value of the spectral index has negligible effect 
on the upper limits. We take into account systematic errors due to uncertainties of the extension, 
the Galactic diffuse model, the LAT effective area, and the effect of the nearby gamma-ray pulsar 
PSR J 1809-2332. Systematic errors associated with the extension are estimated by varying the 
source size by =b 1 a. We also evaluated the effect of the shape by comparing the disk shape 
used in this paper and a circular two-dimensional Gaussian. Systematic errors due to the Galactic 
diffuse model are estimated by using the residual gamma-ray data with respect to the best fit model 
in the region where no LAT source is present, specifically, / = 2.°3-4.°7 andZ7 = -l.°25-0.°75. The 
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observed residual is energy dependent and can be modeled as ~ (116 {El\ GeV)^"*^^'" - 100) % 
of the total Galactic diffuse flux. The normalization of the Galactic diffuse model is adjusted 
according to the above equation to estimate the systematic error on the source flux. We evaluated 
the position dependence of this residual and found a dispersion of ~5%, which is not a large effect. 
Uncertainties of the Galactic diffuse model make a dominant contribution to systematic errors in 
all the energy bins. Uncertainties of the LAT effective area is 10% at 100 MeV , decreasing to 5% 



at 500 MeV, and increasing to 20% at 10 GeV and above (IRando et al. Il2009|) . The effect of the 



nearby gamma-ray pulsar was estimated by varying the flux by ± 1 a. 

We evaluated a possibility of a spectral break in the LAT energy band by comparing the 
likelihood of the spectral fit for the LAT data between a simple power law and a broken power law 
as a spectral model of Source N and Source S. The fits yield the likelihood ratio -21n(LpL/LBPL) ~ 
62 for Source N, where LpL and Lbpl are the likelihoods for the simple power-law model and 
the broken power-law model, respectively. The likelihood ratio slightly decreases to at least ~ 
48 (corresponding to 6.6 a with two degrees of freedom) in the worst case accounting for 1 a 
systematic uncertainties. Thus we conclude that Source N has a spectral break at 1.0 ± 0.2 GeV, 
where the error is dominated by statistics. Photon indices are 2.09 ± 0.08 (stat) ± 0.28 (sys) below 
the break and 2.74 ± 0.06 (stat) ± 0.09 (sys) above the break. Note t hat the photon index a bove 



the break is consistent with the H.E.S.S. measurement of 2.66 ±0.27 (|Aharonian et al.ll2008[) . On 
the other hand, we do not find any evidence of a spectral break for Source S. The photon index is 
found to be 2.19 ± 0.14 (stat) ± 0.41 (sys). 

We placed upper limits on the gamma-ray flux in the LAT band at the positions of HESS J 1 800-240A 
and HESS J1800-240C. The upper limits are compared with the H.E.S.S. spectra in Fig.|4](a) and 
(b). The gamma-ray upper limits in the GeV band at the location of HESS J1800-240A appear 
to be inconsistent if the source spectrum is a simple power law. We evaluated the possibility of 
a break between the LAT and H.E.S.S. energy band by comparing the likelihood of the spectral 
fit for the LAT data between a broken power law and a simple power law as a spectral model of 
this source. We fixed the spectral index and the flux values using the H.E.S.S. measurements for 
the simple power-law model while we fixed the spectral index above the break and the flux at the 
same value as the simple power law and varied the break energy and the spectral index below the 
break for the broken power-law model. We used the parameter values at 1 a away from the best 
fit values from H.E.S.S. in the direction where we expect less flux in the LAT energy band, i.e. 
harder spectral index and lower flux than the best fit values. The fits yield the likelihood ratio 
-21n(LpL/LBPL) ~ 35 (corresponding to a 5.5 a significance with two degrees of freedom) for a 
break energy between 10 GeV and 100 GeV. The likelihood ratio decreases to ~ 8 correspond- 
ing to a ~ 2.4 cr significance if we use parameter values 2 a away from the best-fit values from 
H.E.S.S. We conclude that HESS J1800-240A is likely to have a spectral break between 10 GeV 
and 100 GeV. 
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4. DISCUSSION 

4.1. IFGL J1801.3-2322C (Source N) 

IFGL J1801.3-2322c (Source N) is found to be extended and positionally coincident with 
HESS J1801-233. Given the clear spatial match between the TeV source and the molecular cloud 
interacting with SNR W28 and a smoothly connected spectrum in the GeV-TeV band, we assume 
that the bulk of GeV and TeV gamma-ray emission comes from the molecular cloud illuminated by 
particles accelerated in supernova shocks. Below, we adopt the simplest assumption, that GeV and 
TeV gamma rays are emitted by a population of accelerated protons and electrons distributed in 
the same region characterized by constant density and magnetic field strength. It should be noted 
that the imaging resolution of the LAT is not good enough to differentiate the GeV emission region 
from the TeV region. We assume the injected electrons have the same momentum distribution 
as the protons. This assumption implies a break in the particle momentum spectrum because 
the spectral index of the radio data, corresponding to lower particle momenta, is much harder 
than that for the gamma-rays, which correspond to higher particle momenta. Therefore, we use a 
broken power-law to describe the particle momentum distribution in the region. Electrons suffer 
energy losses due to ionization (or Coulomb scattering), bremsstrahlung, synchrotron processes, 
and inverse Compton scattering. The modification of t he electron spectral distribution due to such 



losses was calculated according to lAtoyan et al.l (| 19951) . where electrons are assumed to be injected 



at t = from an impulsive source. Since diffusive shock acceleration theory generally predicts 
particle accelerations in the Sedov phase with a typical duration 10^-10"* yrs, the assumption of an 
impulsive source would be a good approximation for SNRs with the age of lO'^-lO^ yrs. We adopt 
4 X 10^ yrs for the age of W28 in this modeling. Note that here we consider the minimum momenta 
of protons and electrons to be 100 MeV c"^ since the details of the proton/electron injection process 
are poorly known. The gamma-ray spectrum from tt" decay produced by the interaction of protons 
with ambient hydrogen is scaled by a facto r of 1.84 to account for helium and heavy nuclei in 



target material and cosmic-ray composition (|Morill2009|) . 



First, we consider a rr^-decay model to account for the broadband gamma-ray spectrum. In- 
deed, the hadronic scenario gives the most satisfactory explanation f or the GeV gamrn a rays ob- 



served in other SN Rs interacting with molecular gas such as W51C (lAbdo et al.ll2009|) and W44 



(|Abdo et al.ll2010a|) . The number index of protons in the high-energy regime is constrained to be 
Oh ~ 2.7 from the gamma-ray spectral slope. The observed gamma-ray luminosity requires the 
gas density to be much larger than hh ~ 1 cm"^ averaged over the entire SNR shell in order not 
to exceed the typical kinetic energy of a supernova explosion (~ 10^' erg). The resulting total 
proton energy, Wp ~ 1.3 x 10"^^ • (10^ cnT^ /nu) ■ {d/2 kpc)^ erg, is less than 10% of the typical 
kinetic energy of supernova explosions and quite reasonable. Note that Wp is not the total energy 
of accelerated protons but that of the in-situ protons in the molecular clouds. Using the parame- 
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ters summarized in Table [H we calculated radiation model curves as shown in Figure [5] (a). As 
described above, the spectral index of the proton momentum below the break is determined to be 
Ol ~ 1.7 by modeling the radio spectrum as synchrotron radiation by relativistic electrons because 
we assume that protons and electrons have identical injection spectra. It is difficult to derive the 
break point of the proton momentum spectrum from the break of the gamma-ray spectrum since it 
lies in the region where we expect a gamma-ray spectral curvature due to kinematics of tt" produc- 
tion and decays. Because of this, the gamma-ray spectrum gives an upper bound for the momentum 
break at ~ 5 GeV c"^ . The momentum break cannot be lower than ~ 2 GeV c"^ to avoid conflict 
with the radio data. Here we adopt 2 GeV c~^ . The magnetic field strength is constrained to be 
B ~ 160 fiG, for an electron-to-proton ratio of K^p = 0.01, which is the ratio found in the local 
cosmic-ray abundance. Here the ratio is defined at a particle momentum of 1 GeV c"^ . This choice 
of Kep results in minor contribution from the electron bremsstrahlung emission in the gamma-ray 
band. 

On the other hand, it is difficult to model the GeV-TeV spectrum by the electron bremsstrahlung 
component alone because the break in the electron spectrum corresponding to the gamma-ray spec- 
trum will appear in the radio data as sh own in Fig. [5](b) although it might be contaminated by a 
thermal component (IDubner et al.l 120001) . Moreover, the magnetic field strength is constrained to 
be less than 5^4 /iG for the age of 4 x lO'* yrs since no apparent synchrotron cooling is observed 
at TeV energies. The low magnetic field requires the gas density to be smaller than hh ~ 5 cm"^ 
from the flux ratio of the synchrotron component in the radio band and the bremsstrahlung in the 
gamma-ray band. These low values for the magnetic field and the gas density are in disagreement 
with our assumption that the gamma-ray emission comes from the dense molecular cloud. 

The gamma-ray spectrum is formally reproduced by inverse Compton (IC) emission from 
accelerated electrons in Figl5](c). The interstellar radiation field for the IC scattering (see Table [B 
is comprised of infrared, optical and the cosmic microwave background (CMB). The infrared and 
optical componen ts are taken from the interstellar radiation field at the location of W28 in the 
GALPROP code ([Porter et al.l 120081) . Since the spectral shape of the non-CMB component is 
very complex, it is approximated by two infrared and two optical blackbody components. The 
flux ratio between the IC and the synchrotron components constrains the magnetic field to be less 
than 0.6/iG. The total energy in electrons is calculated to be W^ ~ 9.0 x 10^° ■ (d/2 kpc)^ erg for 
an energy density of ~ 1.8 eV cm"^ for the interstellar radiation field. This model requires a 
rather low gas density of hh ~ 2 x 10"^ cm"^ to suppress the electron bremsstrahlung, which is in 
contradiction with our assumption. 

Observations with Fermi have demonstrated that bright, extended gamma-ray sources co 



incident with middle-aged SNRs interacting w ith molecular clou ds, such as W44 (|Abdo et al 



2010ah . W51C (lAbdoetal.ll2009h . and IC 443 (lAbdo et all 1201 Obi) , exhibit spectral breaks from 
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9 GeV c"^ to a few tens of GeV c"' in the proton momentum spectrum. Our observations of 
IFGL J1801.3-2322c in the vicinity of W28 combined with the radio data constrain the proton 
momentum break to be in the range, 2-5 GeV c"^ . The observed energy distribution of relativistic 
particles (most likely protons) could be influenced greatly by diffusive transport of particles. If 
so, the relationship between the observed particle spectrum and the acceleration spectrum should 
be rather complicated. The steep particle spectrum, a ~ 2.7, deduced for IFGL J 1801.3-2322c 
does not necessarily represent the acceleration index. iGabici & AharonianI (|2007h discussed the 
time evolution of non-thermal emission from molecular clouds illuminated by cosmic rays from a 
nearby SNR and predicted a steep gamma-ray spectrum for an old SNR due to energy-dependent 
diffusion of cosmic rays. More detailed studies of the properties of interactions between SNRs 
and molecular clouds are required for a comprehensive description of the differences in the break 
momenta and photon indices above the break among the above SNRs. 



4.2. IFGL J1800.5-2359C (Source S) 



IFGL J1800.5-2359c (Source S) was found to have no significant extension and spatially 
coincides with the TeV source HESS J1800-240B, which is associated with molecular clouds 
that contain the ultra-compact H II region W28A2 (G5. 89-0. 39). The observed GeV-TeV spec- 
trum for Source S can be formally described by 7r° dominated, bremsstrahlung dominated, and 
IC dominated mode ls in a similar manners to Source N, as shown in Fig. [6l Free-free emission 
(|G6mez et al.lll99l|) is responsible for the radio spectrum of W28A2, which provides only upper 
limits in the SED. Thus v alues of ^ and hh are le ss constrained by pion and bremsstrahlung models. 



We adopt hh = 10^ cm ^ (lAharonian et al. 



2008|) for both models. From the positional coincidence. 



W28A2 can be invoked as a possible source of high-energy particles resp onsible for the gamma-ray 



emission. W28A2 exhibits very ener getic bipolar molecular outflows (IHarvey & Forceillelll988 



Acord et al. 1997: SoUins et al. 



2004|) which would arise from the accretion of matter by a stellar 



progenitor. iKlaassen et al.l (|2006|) estimated the total kinetic energy of the outflow of W28A2 to 
be 3.5 X 10''^^ erg, which requires very high matter density to account for the observed gamma-ray 
flux in both pion and bremsstrahlu ng models. If the matter density is as high as ~ 10^ cm"^ as 
suggested by iKlaassen et al.l (|2006() . the resulting total energy would be less than a few percent of 
the total kinematic energy of the outflow of W28A2. 

Cosmic rays that escaped from the SNR W28 in earlier epochs would be another possi- 
ble explanation given the energetics. Note that Source S would lie at a projected distance of 
~ 10 • (J/2kpc) pc from the southern circular boundary of W28. If we assume that cosmic rays 
are uniformly radiated by the SNR, the cosmic-ray flux bombarding the molecular cloud from the 
SNR can be scaled by the fraction of the solid angle of the molecular cloud, ~ 8 x 10"^. Using this 
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scaling factor, the total energies in the whole SNR W28 could be calculated from the total energy 
observed in this source to be ~ 2 x 10"^^ erg for protons in the tt" model. This is a reasonable 
value for an SNR. The energy density of such cosmic rays is enhanced to be at least ~ 2 eV cm"^ 
at the molecular clouds under the assumption that all of the associated clouds (M ~ 4 x 10"^ Mq 
at 2 kpc) are interacting with energetic particles. The molecular clouds can accumulate those 
cosmic rays in the past since the diffusion coefficient is expected to be low in dense environ- 



ments ( Aharonian & Atovanlll996 ). 



While we detected LAT sources spatially coincident with HESS J 1801-233 and HESS J1800-230B, 
no strong LAT counterpart can be found at the positions of HESS J1800-230A and HESS J1800-230C. 
We found evidence of a spectral break between 10 GeV and 100 GeV for HESS J1800-230A. The 
velocity with respect to the local standard of rest of the molecular emission po sitionally coinci- 



dent w ith this source peaks at ~16 km s \ corresponding to a distance of ~4 kpc (|Aharonian et al 



20081) . However, the velocity difference between this cloud and the one associated with HESS J 1 800-230B 
is only ~ 10 km s"^ and can be attributed to the proper motion of the cloud complex that contains 
both clouds. If the distance to the cloud associated with HESS J1800-230A is 4 kpc, then it is 
unlikely that this source is associated with W28 since the distance to W28 is ~2 kpc. H U regions, 
G6. 1-0.6 and/or 06. 225-0. 569 could be a source of accelerated particles. If the cloud associ- 
ated with HESS J1800-230A is part of the cloud complex that includes the cloud associated with 
HESS J1800-230B, then the distance would be 2 kpc, which is very similar to the clouds associ- 
ated with W28. In this case, particles escaped from W28 may be responsible for the gamma-ray 
emissions at both HESS J1800-230A and B and further studies of the origin of the difference in 
the spectral shapes in these sources may provide good constraints on the particle diffusion process. 



5. CONCLUSIONS 

We have investigated two LAT sources in the W28 field. Source N (IFGL J1801.3-2322c) 
which is located at the northeast boundary of the SNR W28 is positionally coincident with shocked 
molecular clouds, and is spatially extended. The spectrum has a break around 1.0 GeV and 
smoothly connects to the TeV spectrum, suggesting a physical relationship. Decay of 7r°s pro- 
duced by the interaction of an SNR with molecular clouds naturally explains the gamma rays from 
Source N based on the spatial correlation between GeV gamma rays and molecular clouds and 
the energetics of cosmic rays. Electron bremsstrahlung can not be ruled out completely although 
it requires a low density and low magnetic field in contradiction with the association with the 
molecular clouds. W28 is the most plausible energy source due to the observational evidence of 
interaction with the molecular clouds. The soft spectrum of the gamma rays may be explained 
by the time evolution of non-thermal emission from molecular clouds illuminated by cosmic rays 
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from a nearby SNR due to energy-dependent diffusion of cosmic rays. 

The Source S (IFGL J1800.5-2359c) was found to have no significant extension and spatially 
coincides with the TeV source HESS J1800-240B, molecular clouds, and the ultracompact H II 
region W28A2 (G5.89-0.39). The compact H II region W28A2 and the SNR W28 are possible 
energy sources, but the W28A2 hypothesis requires extremely dense gas. 

While no significant LAT counterpart is found at the positions of HESS J1800-230A and 
HESS J1800-230C, the LAT upper limits for HESS J1800-230A coupled with the H.E.S.S. data 
points imply a spectral break between 10 GeV and 100 GeV. 
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Fig. 1 . — Fermi LAT 2-10 GeV counts map around the supernova remnant (SNR) W28. The count 
map is smoothed by a Gaussian kernel of a = 0.°2, with the pixel size of 0.°025. The white line 
from top left to bottom right indicates the Galactic plane. 
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Fig. 2. — Close-up views of the LAT 2-10 GeV count map around W28. The counts map is 
smoothed by a Gaussian kernel of cr = 0.°2, with the pixel size of 0.°025. The inset of each figure 
shows the effective LAT PSF for a photon spectral index of 2.5. A green circle in the north of 
each figure indicates the best-fit disk size for Source N. A green cross indicate s the position of 
Sourc e S. White diam onds indicate H II regions; W28A2 (see text), 06. 1-0.6 (IKuchar & Clark 
1997k G6.225-0.569 (lLockmanlll989l) . The diamond on the right is W28A2. Black contours in 
(a) sh ow the H.E.S.S . significance map for TeV gamma rays at 20, 40, 60 and 80% of the peak 
value (JAharonian et a l. 2008). Bright TeV spots in the south are HESS J 1800-240 A, B and C as 
indicated in the figure. Black contours in (b) give CO (J=l-0) line intensity taken by NANTEN at 
25, 50, 75 % levels, for the velocity range f rom km s~^ to 20 km s~\ corresponding to k inematic 



distances of approximately to ~ 4 kpc (JMizuno & Fukuil 120041: iTakeuchi et al. 



2010). Black 



contours in (c) indicate the VLA 90cm image at 25, 50, 75 % of the peak intensity (JBrogan Il2006|) 
Outer boundaries of SNRs, as determined by the radio images, are drawn as white dashed circles 
A white plus sign shows the position of PSR J 180 1-23. 
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Fig. 3. — (a) Spectral energy distribution of the Fermi LAT source on the northeast boundary of 
the supernova remnant W28. The red squares in the GeV regime are the LAT data. Horizontal 
bars for data points indicate the energy range used in the fit. Vertical bars show statistical errors 
in red and systematic errors in black. Upper limits are obtained at the 90% confidence level in 
energy bins where the likelihood ratio va lue is < 9. The black circles represent data points for 
HESS J1801-233 jAharonian et allboosl) . (b) Spectral energy distribution of the LAT source to 
the south of W28. The red region is the 68% co nfidence range of the LAT spectrum. The black 
circles show data points for HESS J1800-240B (lAharonian et al.ll2008k 
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Fig. 4. — Upper limits obtained from the LAT data at the 90% confidence level in each energy 
bin at the positions of (a) HESS J1800-240A and (b) HESS J1800-240C, on the assumption of 
a phot on index of 2 for the power-law function. Black circles are H.E.S.S. spectral measure- 
ments ( Aharonian et al.ll2008b . 
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Fig. 5. — Multi-band spectra of the Fermi LAT source on the northeast boundary of the supernova 
remnant W28. The red squares in the GeV regime are the LAT data, where t he red and black 



errors on the flux ar e statistical and systematic, respectively. The radio emission (iKovalenko et al 



19941 : iDubner et al.l 12000 ) is modeled by synchrotron radiation, while the gamma-ray emission 
is modeled by different combinations of 7r°-decay (long-dashed curve), bremsstrahlung (dashed 
curve), and inverse Compton (IC) scattering (dotted curve). Details of the models are described in 
the text. 
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Fig. 6. — Multi-band spectra of the Fermi LAT source to the south of the supernova remnant W28. 
The red squares in the GeV regime are the LAT data, where the red and black errors on the flux 



are st atistical and systematic, respectively. The radio emission (JHanda et al.lll987l : iGomez et al 



19911) is modeled by synchrotron radiation, while the gamma-ray emission is modeled by differ- 
ent combinations of 7r°-decay (long-dashed curve), brems strahlung (dashed curve), and inverse 
Compton (IC) scattering (dotted curve). The sum of the three components is shown as a solid 
curve. Details of the models are described in the text. 
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Table 1: Model parameters for the Fermi LAT sources around W28. 



Model 


Kep^ 
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Pb^ 
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(GeV c-i) 




(fiG) 


(cm-3) 


(10^" erg) 


(lO^'^erg) 


Source N 


















(a) Pion 


0.01 


1.7 
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2.7 
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100 


1.3 
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(b) Bremsstrahlung 
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The ratio of electron and proton distribution functions at 1 GeV c ' . 

''A momentum distribution of particles is assumed to be a broken power-law, where the indices and the break momentum 

are identical for both accelerated protons and electrons. The q;l is a photon index below the momentum break. 
"^The pb is a momentum break for particle distribution. 

''The photon index for broken power-law functions above the momentum break. 
'"Average hydrogen number density of ambient medium. 

^The distance is assumed to be 2 kpc. The total energy is calculated for particles > 100 MeV c"'. 
*Seed photons for inverse Compton scattering of electrons include the CMB, two infrared (Tjr = 29,4.9 x 10^ K, 

t/iR = 0.29,5.3 X 10"^ eV cm"^, respectively), and two optical components (Topt = 3.6 x 10-^, 1.0 x 10"* K, f/opt = 

0.37,0.13 eV cm"-', respectively) in the vicinity of W28, assuming a distance of 2 kpc. 



